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Abstract 

We consider radiative and pionic decays of the new D s -mesons in the framework of 
a phenomenologically motivated approach. Present data on ratios of the two kinds of 
decays can be described without explicit using a 4-quark component. Most probably, 
the isospin violation in decays of different D s -mesons is not universal, and the binding 
potential should be different from Coulombic. New precise measurements may provide 
further clarification for the nature of the excited states. 



1 Introduction 



Recent discoveries of several new charmed meson states are among the most impressive 
recent events in the heavy flavour physics. These are, first of all, two narrow charmed- 
strange states, D s j(2317) and _D s j(2460) . The former one was initially observed by the 
BaBar group PQ, and later confirmed by CLEO J2j and Belle jHIHj- Some evidence for the 
second state was also noted in ref. [T], but it was seen definitely only later, by CLEO [2J 
and Belle 011]. Further, it was reliably confirmed by BaBar as well [5 j. To the present 
moment the two states have been seen both in B-decays and in the continuum annihilation 
e + e _ — > cc j3J. So they are definitely established, though their properties, especially their 
quantum numbers, are not quite clear yet (experimental minireviews on these states see 
in 0[7j). 

Two more states, non-strange charmed mesons, seem to be discovered by Belle |H] in 
final states D^n^ and D* ± 7r T . They have been extracted from the coherent amplitude 
analysis of the Dalitz plots for Dim and D*7itt . As a result, their quantum numbers 
are fixed (dominant S- waves), but their existence strongly needs confirmation. It is even 
more so, that these states have very large widths (some hundreds MeV), and overlap 
with kinematical reflections from other resonances. Evidence for the broad 5*-wave state 
in the 2?7r-system has been recently presented also by the Wideband photoproduction 
experiment FOCUS [HJ. However, the masses as given by the two collaborations differ by 
a couple of standard deviations, and the situation needs further clarification. 

Unexpected properties of the new Instates have induced active theoretical discus- 
sions. Presumably, the mesons have spin-parity numbers + for D s j(2317) and 1 + for 
D s j(2A60) (these values, at least, do not contradict any experimental data). States with 
such quantum numbers have been, of course, predicted earlier, but they were expected 
to have higher masses and much higher widths |1L)| II lj. Therefore, various suggestions 
have appeared in the literature that the new states may have some nature different from 
the canonical qq states. The published ideas include Di^-molecules [12} . Dn-atom [T3], 
4-quark states [11} . and their combination with qq states 15 . 

In the framework of more "usual" approaches, the chiral perturbation theory was 
combined with HQET to include the new states into D s spectroscopy as chiral partners 
of the well-known states D s (0~) and D*(l~) [TH } I17 [ IT%}. The previous calculations for 
masses were mainly based on the relativistic quark model with the potential consisting 
of two parts: Coulomb interaction at small distances (due to one-gluon exchange), and 
linear potential at large distances (for confinement) 1201 E] Now it was suggested 
that such method could be applied to masses of the new states as well, if one modified 
the confinement potential [22J- 

Decay properties of the new states have been also considered in various assumptions, 
e.g., HQET with chiral effective Lagrangian for light quarks fniEBJ, previous |24J and 
modified [22] potential pictures. However, all they show disagreement with present exper- 
imental data, especially for the ratio of radiative and hadronic decays of the new mesons. 
This has led to the rather pessimistic conclusion that if the data are confirmed "the 
conventional D* and Dl{ 2 states have yet to be discovered" [2~lj . 
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In this note we use another, essentially phenomenological way to estimate partial 
widths of the new D s states. We try to use as few dynamical assumptions as possible. 
Our approach is explained in Section 2. Then it is applied to the new meson decays 
(Sections 3,4). Our results are summarized and discussed in Section 5. 

In what follows we always assume that the states D s j(2317) and _D s j(2460) have indeed 
the spin-parity + and 1 + respectively. We will denote them as _D* (2317) and _D s i(2460) . 



2 Mass dependence of matrix elements 

To explain our approach, let us discuss the problem in terms of 2-body configurations 
with an efficient potential. Then decay amplitudes are expressible through some overlap 
integrals between initial and final wave functions. 

For the pure Coulomb potential, having no specific mass-dimensional parameter, the 
overlap integrals could be written as a power of the reduced mass \i = m 1 vn%j (m 1 + m 2 ) 
multiplied by a universal dimensionless factor (independent of the constituent masses). 
The same would be true, in particular, for mass differences of radial excitations having 
the same flavour quantum numbers. Those mass differences would be proportional to the 
corresponding reduced masses. 

Experimental data, however, do not support such expectations. Let us take, for in- 
stance, vector and tensor mesons with different flavour content. They are well established 
and well studied for all known flavours (excluding beauty mesons), they nicely correspond 
to expectations for the quark- ant iquark system. The both states have the same internal 
spin structure (the total spin of constituents is S = 1), and differ only in the internal 
orbital momentum (L — 0, 1 for the vector and tensor states respectively). Moreover, 
tensor and vector non-strange isospin-singlet mesons have nearly the same SU(3)p vio- 
lating mixing angles, leading to nearly ideal separation of the strange and non-strange 
components. It is reasonable, therefore, to compare tensor-vector mass differences for 
mesons of various flavours. 



For the mesons with light (u,d,s) quarks we have the differences 

M a2 -M p = (547 ± 1) MeV , (1) 

M K * - M K * = (534 ± 2) MeV , (2) 

M h -M u = (493 ± 1) MeV , (3) 

Mf 2 -M+= (505 ± 5) MeV , (4) 

which are nearly the same. Note that the change (u, d) — > s , though essentially shifts the 
particle masses, produces only very small effect (~ 2.5%) onT-y mass differences. 

For mesons with heavy quarks the differences are j25j 

M D ,o - M D *a = (452 ± 2) MeV , (5) 
M D *± - Md»± = (449 ± 4) MeV , (6) 
M D *± - M D ,± = (460 ± 2) MeV , (7) 
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M Xc2 - Mj/tp = (459.31 ± 0.14) MeV , (8) 
M Xb2{1P) - M T{15) = (452.3 ± 0.6) MeV . (9) 

Again, the change of the mass difference under substitution (u, d) — > s is only ~ 2.5% . 
With this accuracy the T — V mass differences are the same for all mesons containing 
at least one heavy quark, c or b, being lower (!) than ones for light-quark mesons 1 . It is 
interesting to note that in the sequence of hidden- flavour quarkonia [/2,k>] — > [f^fi] ~^ 
[Xc2,J/ip] — ► [Xb2{^-P)i T(IS')] , which corresponds to the quark sequence u/d — > s — > 
c —>■ b with monotonously growing constituent masses, the T — V mass difference slightly 
increases at the first step, and then monotonously decreases. Such behaviour does not 
correspond to Coulomb regularities, even for the most compact, 6-quark systems. 

The discussed T — V mass differences are free from the constituent "rest mass" con- 
tributions, and are expressed only through diagonal matrix elements of some operators 
(corresponding to the kinetic and potential energy) . As we see, such matrix elements, for 
changing quark contents, vary not more than ~ 20% , i.e., they reveal only low dependence 
on the quark masses. 

Decay amplitudes are also expressed through some matrix elements, though nondiago- 
nal ones, and of different operators. Nevertheless, we will assume that the emerging matrix 
elements also have low dependence on the quark masses. As the first approximation, we 
will neglect such dependence. 



3 Decays of L> s0 (2317) 

Let us begin with consideration of _D* (2317) . Its only observed decay mode is 

D£(2317) ^Dfn\ (10) 

having the final momentum about 320 MeV/c. If the initial meson corresponds to the 
cs system, then this decay violates isospin, and should have some suppression. It is well 
known that there are at least two possible sources for such violation, electromagnetic 
interactions and/or mass difference of d- and li-quarks. Examples of the isospin violation 
for light-quark hadrons (without explicit photon radiation or absorption) demonstrate 
numerical smallness of order 10~ 2 — 10~ 3 ~ 0(a) (the rare exclusion, rather intensive 
decay to — > 7r + 7r~ , only confirms this empirical rule, since contains enhancement due to 
the small puj mass difference). This does not allow to separate strong and electromagnetic 
violations of the isospin. 

Situation for heavy-quark hadrons seems to be different. The best studied example is 
the isospin-forbidden decay ip{2S) — > J/ipn with branching ratio 

B n o = (9.6 ±2.1) ■ 10~ 4 

1 Recall that T- and T^-mesons differ only in the internal orbital momentum. Mass difference for such 
systems bound by the Coulomb potential should grow with increase of constituent mass(es). In contrast, 
V- and P-mesons differ only in spin orientation, and their mass differences should decrease with increas- 
ing constituent masses, independently of the binding potential, due to decrease of (chromo)magnetic 
moments. This latter expectation exactly agrees with measurements. 
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and final momentum p w o = 527 MeV/c. It can be compared with the similar, but isospin- 
allowed decay ip(2S) — > J/iprj having branching ratio 25 

B v = (3.17 ±0.21) • 1(T 2 

and the final momentum p v = 200 MeV/c. Both decays produce p-wave final states. So, 
the suppression factor in the isospin-violating amplitude is 

1/2 

= (4.07 ±0.47) ■ 10~ 2 . (11) 

Such a value looks too large for the electromagnetic mechanism. It might be enhanced by 
7] 7r°-mixing, in similarity (though weaker) with p° cu-mixing which enhances the isospin- 
violating decay uj — > 7r + 7r~. If so, the same parameter e should work in all decays which 
produce ir° with violation of isospin. Alternatively, its large value in compact charmonium 
mesons could give evidence for the direct effect of the mass difference — m u , being 
determined, say, by the ratio (m^ — m u ) j Aqcd • In any case, we will try to assume that 
this parameter is universally applicable to isospin violation in decays of both charmonium 
and Z} s -mesons. 

Now we can compare decay (jlOj) with the similar, but isospin-conserving decay 

D*± -> D ± 7T° , (12) 

where D^ 1 are charmed-nonstrange mesons with J p = + . Their isotopic partner, Dq° , 
was found by Belle [H] to have the mass and total width 

M* = (2308 ± 37) MeV , T* = (276 ± 66) MeV . (13) 

With reasonable accuracy, we can assume for the same values of the mass and total 
width. Further, we suggest that the large total width is completely due to the decay mode 
Dq — ► Dtt which should be strongly dominant. For the particular charged channel (112)) 
we then obtain the final momentum about 378 MeV/c and the partial width 

T(L)* ± -> D ± 7T°) w (l/3)r* = (92 ± 22) MeV . (14) 

The rest (2/3)Tq comes, according to isotopic relations, from the mode Dq — > D° -k^ . 

The both decays ()10|) and (|12|) have S"-wave final states; each of their widths should 
be a product of the final momentum by the square of some matrix element. We can 
assume that those matrix elements differ only by the suppression factor e of Eq. (jllj) . 
Then, accounting for this suppression and for the increased final momentum, 378 MeV/c 
in the decay (fT2"j) vs. 320 MeV/c in the decay (jTDjl. we estimate, at last, 

r(D% -> Df vr°) = (129 ± 43) keV . (15) 

On the other hand, according to FOCUS the enhancements near the mass 2300 MeV 
are consistent with kinematical effects of feed-downs from different resonant states, while 
the true scalar states Dq° and Dq have 

M* = (2407 ± 41) MeV , T* = (240 ± 81) MeV ; (16) 
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M* ± = (2403 ± 38) MeV , T* ± = (283 ± 42) MeV . 



(17) 



At such parameters the decay (|T2*|) has the final momentum 459 MeV and leads to the 
estimate 

T(D*± -> Df vr°) = (109 ± 16) keV . (18) 

It is somewhat lower but consistent with the Belle-based estimate (|15|). 

Let us compare now the pion decay (fTUj) with radiative decays of Dfj(23l7) . None of 
such radiative decays has been observed yet. If the Dfj(2317) indeed has J p = + , then 
its decay to D s , with J p = 0~ , is strictly forbidden (as any radiative 0-0 transition). 
Possible is the decay 

1^(2317)^/^7. (19) 

Its absence in experiment is one of strange properties of -D s j(2317) . This decay should 
exist, if the interpretation of D s j(2317) as D* is correct. To estimate its expected prob- 
ability, we compare it to the similar decay 

X c o->W7- (20) 

The both decays should correspond to El transitions. Their partial widths may be 
written as 

r(o + ^i- + 7 ) = |(d)| 2 ^. (21) 

Here (d) may be interpreted as the transition matrix element of the operator of electric 
dipole moment, p 7 is the momentum of the produced photon. In its turn, the dipole 
moment element may be written as 

(d) = (Ae) ■ (r) , (22) 

where (Ae) is the weighted charge difference of constituents in the 2-body (qq) system. 
For the charmonium (cc) it is 

(Ae) CJ = (e c - es)/2 = 2/3 . 

For the -D s -system (cs) , with unequal quark masses, we take non-relativistic expression 

(Ae) C s = (m c &r - m s e c )/{m c + m s ) . 

Using the naive estimation m c « 3m s , we obtain (Ae) C s ~ (Ae) C c/8 . Fitting experimental 
meson spectra to some detailed potential calculations gives m c w 4m s (see, e.g., ref . |2~4] ) : 
then (Ae)cs ~ (Ae) C c/5 . So we can safely take 

(Ae) C s<l(Ae) c - c . (23) 
5 
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The matrix element (r) could also be different for the systems cc and cs (in the pure 
Coulomb case the difference would be determined by the factor (m c /m s + l)/2, which 
can be ~ 2 for realistic constituent masses). However, even ground state of charmonium 
is not sufficiently compact to be concentrated in the pure Coulomb area. As a result, 
this matrix element should change much smaller. For simplicity, as in other cases, we 
take it independent of the constituent masses, and apply relation (J2~3j) directly to matrix 
elements of the dipole moment. 

The partial width of decay (j20j) is known experimentally to be [23] 

r(Xco -> J/tpj) = 119 ± 19 keV. (24) 

Its final photon momentum is 303 Mev/c, vs. about 200 MeV/c in decay (|19|). Accounting 
for suppression (|2*3*|) and for the p 3 behaviour of Eq. ()21)1 with the smaller final momentum, 
we can now estimate expected partial width of the decay 

r(D* s0 ^D* sl )<lAkeV. (25) 

Treatment of Eq.(|23j) as equality with the "naive" coefficient 1/8 instead of "safe" 1/5 
would give the value 

r( J D: ^ J D: 7 )^0.55keV. (26) 



Thus we expect the ratio T(D* — ► D*^f)/T(D* — > D s tt°) to be not higher than 
~ 1.5% . The present experimental boundary of CLEO jBj 

13^(2317)-/^°) < 6% (27) 
completely agrees with our expectation. 



4 Decays of L> s i(2460) 

Here we begin with considering the radiative decay 

£> sl (2460)^£> s7 , (28) 

which was definitely observed by Belle jlj. The Z) s -meson has J p = 0~, and for D s i(2460) 
we assume J p = 1 + . Then the radiative decay (|2*%|) corresponds to El transition, and can 
be estimated in essentially the same way as done in the preceding section for the radiative 
decay of D* . 

Charmonium analog of D s is r] c with J PC = h . The state -D s i(2460) has two 
possible analogs, h c with J PC = l + ~, and Xci with J PC = 1 ++ . Decay Xd ~ ¥ Vcl is 
forbidden by C-parity, and thus, analogous to decay (J28j) could be only the radiative 
decay h c — > rjc'j . However, the state h c itself is badly known (see Listings in tables |2*5]). 
and its radiative decay has never been observed. Instead, we again can use decay P0|. but 
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with an essential note. Contrary to decays (j!9|) and (|20|). which are + — > 1 transitions, 
decay (}2"%j) corresponds to the transition 1 + — > 0~. This leads to the radiative width 

r(l + -0-+ 7 ) = i|(d)| 2 4, (29) 

which should be compared with the expression (}2~T]) . The new factor 1/3 appears due to 
averaging over spin 1 states of the initial system, instead of summing over them for the 
final system in expression ()21|). Taking the other factors to be the same as before, and 
accounting for the larger final momentum in decay (|2"K|) . about 440 MeV/c, we obtain the 
upper boundary 

T(D sl ^D sl ) <5keV. (30) 

It would give 

r(D sl ^D sl )&2keV (31) 

at the "naive" treatment of Eq. ()23)1 as an equality with the coefficient 1/8, instead of 1/5. 
Pionic decay 

L> s i(2460) ^D*tt°, (32) 

with = 425 MeV/ c, violates isospin symmetry. Its theoretical consideration may appear 
ambiguous since there are two possible pion transitions 1 + — > 1~ , through S- or ZJ-wave. 
Of course, they can be discriminated experimentally, but this has not been done yet for 
£ s i(2460). 

In the charmed-nonstrange sector there are two l + -states [BHHH^H] with nearby masses. 
One of them is relatively narrow, another is much wider. Those are D\ with 

M Dl = (2421.4 ± 1.7) MeV , T Dl = (23.7 ± 4.8) MeV (33) 

and D[ with 

M D[ = (2427 ± 36) MeV , r D / = (384 ± 120) MeV . (34) 

We use here the two masses and widths from the simultaneous coherent analysis of the 
Z}*7r-system jS]. PDG-tables contain only a narrower state, its parameters being 
consistent with values (J33|) . On the other hand, analysis of the FOCUS Collaboration [0] 
reveals only a wider state, which mass and width are consistent with (J34j) . 

Further, according to Ref.|H], D\ decays mainly through D-w&ve, and D[ through S- 
wave (the latter agrees with data [Oj), thus explaining the large difference of the widths. 
The mutual admixtures are about 10% in the wave functions [H], and as the first step we 
neglect the mixing. 

Here we need to make an assumption about partial- wave properties of the decay (|3*2*|) . 
If it goes (mainly) in S"-wave, we should compare it to the S"-wave decay D[ — > D* tt° , 
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with the final momentum 482 MeV/c. This can be done in exactly the same way as for 
decays (fTUj) and (fT2"j). and leads to 

r (s) (D sl (2460) -> D* 7t°) = (187 ±73) keV . (35) 

But if the decay (|3~2"j) is D-wave one, it should be compared to the D-wave decay Di — > 
D* 7T° , having the final momentum 477 MeV/c. The width for such decays may be written 
as product of the 5-th power of the final momentum by the square of a matrix element. 
As before, we relate the matrix elements for isospin- violating (D s \) and isospin-conserving 
(.Di) decays by the suppression parameter e and then account for the p 5 -behaviour of the 
partial width. Thus we obtain 

T (D) (D sl (2460) -> D* s 7T°) = (7.4 ± 2.3) keV . (36) 



To discriminate between these two possibilities, we can use the experimental fact of 
observing decay (|2*H|) with relative intensity [HI H] 

r(fl fll (2460)->D,7) n , ( , 7] 
r(£> sl (2460) ^,D s *7r ) 1 ' 

Comparison of values (J30|) . (}35j) and (}36j) shows that our estimate for the S- wave case 
leads to strong disagreement with the ratio (p?7|) . while it is quite compatible with the 
.D-wave case. 



5 Discussion and conclusion 

Let us summarize our assumptions and their implications. 

1) First of all, we have assumed that overlap integrals in decay amplitudes are nearly 
independent of the constituent quark masses. This assumption would be incorrect if the 
binding potential were Coulombic. But realistic potential should be more complicated, 
thus making our assumption admissible. In any case, it is inspired by the phenomenology 
of measured meson mass differences. 

2) Another assumption may be less clearly motivated. It is the universal character of 
the suppression factor e, as given by Eq.QUJl. for isospin- violating amplitudes in decays 
of mesons containing at least one heavy quark. We will discuss this point below in more 
detail. 

All our estimates begin with decays of mesons which are usually considered as quark- 
antiquark systems (charmonium, excited D-mesons). Therefore, we never use, at least 
explicitly, presence of 4-quark components. 

With such assumptions we are able to describe the present data on relation of photon 
and pion decays of the new D s mesons, but only if the decay (J3~2"j) . i.e. -D s i(2460) — > D* 7r°, 
goes (mainly) through D-wave. In terms of the effective heavy quark description, this 
would mean that D s i(2460) corresponds to the state with j = 3/2 . Though the question 



8 



can (and will) be solved by direct measurement of the angular distribution in the decay, 
today such prescription disagrees with familiar expectation of the theoretical community. 
Moreover, if it is correct, another assumed 1 + state, _D s j(2535) , should have j — 1/2 and 
decay to D*K through S'-wave. In that case one could hardly understand its small width 
Y tot < 2.3 MeV |25]. 

If _D s i(2460) corresponds to j = 1/2 , our calculations can still be made consistent with 
data, if we weaken the above assumptions. As the easiest way to this goal we can drop out 
the universality of the suppression factor (fTTj) for isospin violation. Indeed, if we accept 
estimates (j2T?j) and (}3*Tj) for radiative decays, then the amplitude for the S'-wave decay 
(|3^|) should be suppressed by the factor ~ e/7, instead of ~ e , to satisfy experimental 
relation ()37j) . Such suppression (~ 6 • 10~ 3 ) has rather familiar order of smallness and 
does not look too severe. 

On the other side, to satisfy the boundary (|2Tj) . the amplitude for decay ((TOJ) , i-e. 
/}*,^(2317) — > Dfir , should not be suppressed stronger than ~ e/4, which is ~ 10~ 2 . 
Thus, the isospin violation in decays (fTUj) and ()32|) seems to be non-universal, and different 
from violation effects in decays of charmonium. 

This might look strange, since all the cases should have universal contributions due to 
(jl 7r°)-mixing. Note, however, that for charmonium we compare decays to rj and to 7r°, 
while in estimating decays (fTUj) and (|3~2"j) we compare pionic decays of different systems, 
cs and cd . Such ratios should not be necessarily the same. 

Moreover, as shown in Ref.[26j, the universal mixing violation of isospin symmetry 
should be always accompanied by non-universal direct violation (compare the mixing vs. 
direct violation of CP-parity, say, in kaon decays). Due to the small mass difference, the 
(a; p°)-mixing strongly enhances violation, and makes the whole effect to be nearly uni- 
versal. The enhancement due to (rj 7r°)-mixing is weaker, and non-universal contributions 
of direct violation may become more essential. Thus, non-universality of isospin violation 
in decays of different _D s -mesons is not amusing. 

One may also try to weaken the first of our above assumptions. It is easy to see that 
some increase of the overlap integrals for D s -mesons as compared with charmonium allows 
to reach better agreement between our estimates and experimental data. It could be done 
even with universal suppression factor of isospin violation, but then the increase should 
be faster than for the Coulomb potential, in contradiction with data on mass differences. 
Therefore, the suppression, most probably, should be non-universal. We will not discuss 
here structure of the admissible effective potential. 

These two modifications of our assumptions provide interesting features. Estimates 
(J2HJ) and ()31|). appended by non- universal isospin violation, provide for the new D s mesons 
the total widths of less than ~ 10 keV. Universal isospin-violation suppression, together 
with increasing overlap integrals, would give the total widths ~ 100 keV. So, the two 
possibilities could be discriminated by special precise experiments. 

As a conclusion, we can say that present data on decay properties of new D s mesons 
may be described in terms of a quark- ant iquark bound system with some non-Coulombic 
potential. Presence of an essential 4-quark component is not required. Similar conclusion 
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was derived earlier [22] on the basis of spectroscopy only. Our additional, and new, result 
is probable non-universality of isospin violation in pionic decays of different D s states, 
even if it is related to (rj 7r°)-mixing. Important and clarifying new information could 
come from measurements of total widths for the new mesons, though this seems to be a 
very hard problem. 

Acknowledgments 

The work was partly supported by the Russian State Grant SS-1124.2003.2, and by the 
Sofja Kovalevskaja Programme of the Alexander von Humboldt Foundation, the German 
Federal Ministry of Education and Research, and the Programme for Investment in the 
Future of the German Government. 

References 

[I] B.Aubert et al. [BaBar 
|hep-ex/0304020| l]. 

[2] D.B esson et a l. [CLEO 
|hep-ex/0305100| . 

[3] Y.Mi kami et a l. [Belle 
|hep-ex/0307052| . 

[4] P.Krokovny et al. [Belle 
|hep-ex/03080T9| . 

[5] B.Aubert et al. [BaBar Collaboration], submitted to Phys.Rev. D |hep-e x/03 10050 1 . 

[6] F.C.Porter, to appear in Proceedings of the International Europhysics Conference on 
High Energy Physics, July 17-23, 2003, Aachen |hep-ex/ 031 2019] . 

[7] J.Wang, to appear in AIP Proceedings of Beauty 2003, Carnegie-Mellon University, 
Pittsburg, Oct. 14-18, 2003 hep-ex/ 031 2 039| . 

[8] K.Abe et al. [Belle Collaboration], to be submitted to Phys.Rev. D |hep-ex/0307021 . 

[9] J.M.Link et al. [FOCUS Collaboration], submitted to Phys.Lett. jhep-ex/0312060| . 

[10] J.Bartelt and S.Shukla, Ann.Rev.Nucl.Part.Sci 45, 133 (1995). 

[II] M.DiPierro and E.Eichten, Phys.Rev. D64, 114004 (2001). 

[12] T.Barnes, F. E.Close and H.J.Lipkin, Phys.Rev. D68, 054006 (2003) 
|hep-ph/0305025] . 



Collaboration] , 

Collaboration] , 
Collaboration], 
Collaboration], 



Phys.Rev.Lett. 90, 
Phys.Rev. D68, 
Phys.Rev.Lett. 92, 
Phys.Rev.Lett. 91, 



242001 (2003) 

032002 (2003) 

012002 (2004) 

262002 (2003) 



10 



[13 
[14 

[15 
[16 

[it; 
[is; 

[19 
[20 
[21 
[22 
[23 
[24 
[25 



A.P.Szczepaniak, Phys.Lett. B567, 23 (2003). 

H.Y.Cheng and W.S.Hou, Phys.Lett. B566, 193 (2003). 
R.Terasaki, Phys.Rev. D68, 011501 (2003). 



T.Browder, S.Pakvasa and A.Petrov, Phys.Lett. B578, 365 (2004) [hepjph/0307054 . 

W.A.Bardee n, E.J.Eichten and C.T.Hill, Phys.Rev. D68, 054024 (2003) 
|hep-ph/0305"049| 

M.Nowak, M.Rho and I.Zahed, |hep-ph/0307102 



M.Nowak, talk at Hadron'03, Aschaffenburg, Germany, Aug.31-Sep.6, 2003 
|hep-ph/0311154| . 

A.Deandrea, G .Nardulli and A.D.Polosa, Phys.Rev. D68, 097501 (2003) 
|hep-ph/0307069| . 

S.Godfrey and N.Isgur, Phys.Rev. D32, 1889 (1985). 
S.Godfrey and R.Kokoski, Phys.Rev. D43, 1679 (1990). 
E.Eichten, C.Hill and C.Quigg, Phys.Rev.Lett. 71, (1993) 4116. 
R.N.Cahn and J.D.Jackson, Phys.Rev. D68, 037502 (2003) |hep-ph/03050T2| . 
PColangelo and F.DeFazio, Phys.Lett. B570, 180 (2003) |hep-ph/0305l"40l . 
S.Godfrey, Phys.Lett. B568, 254 (2003) [hep-ph/0305T22| . 

K.Hagiwara et al. [Particle Data Gr oup], Phys.Rev. D 66 010001 (2002) and 2003 
partial update for edition 2004 (URL:http://pdg.lbl.gov). 

[26] Ya.I.Azimov, Eur.Phys.J. A16, 209 (2003) |hep-ph/0209T53| . 

Ya.I.Azimov, Ya dern.Fiz. 66, 2039 (2003); Phys.Atom.Nucl. 66, 1989 (2003) 
|hep-ph/0301224| . 



11 



